COMPUTING exp(−τA)b WITH LAGUERRE POLYNOMIALS
نویسندگان
چکیده
Abstract. This paper discusses a method based on Laguerre polynomials combined with a Filtered Conjugate Residual (FCR) framework to compute the product of the exponential of a matrix by a vector. The method implicitly uses an expansion of the exponential function in a series of orthogonal Laguerre polynomials, much like existing methods based on Chebyshev polynomials do. Owing to the fact that orthogonal polynomials satisfy a three-term recurrence, what these series expansion methods offer over other approaches such as Krylov subspace methods lies in the elimination of inner products and the economy in storage since there is no need to compute and keep a set of basis vectors. Compared with Chebyshev polynomials that are orthogonal within a restricted interval and need estimates of the outermost eigenvalues, Laguerre polynomials offer the added feature that they are orthogonal on the half real line, alleviating therefore the need to estimate eigenvalues.
منابع مشابه
Computing exp ( − τA ) b by the Filtered Conjugate Residual Algorithm ∗
This paper discusses a class of Filtered Conjugate Residual Algorithms (FCR) as a way to compute the product of the exponential of a matrix by an arbitrary vector. These methods utilize implicitly expansions of the exponential function in a basis of orthogonal polynomials. FCR methods based on Laguerre, Hermite, and Chebyshev polynomials, are described and their performances are compared. The p...
متن کاملAddition Theorems as Three-Dimensional Taylor Expansions. II. B Functions and Other Exponentially Decaying Functions
Addition theorems can be constructed by doing three-dimensional Taylor expansions according to f(r + r) = exp(r · ∇)f(r). Since, however, one is normally interested in addition theorems of irreducible spherical tensors, the application of the translation operator in its Cartesian form exp(x∂/∂x) exp(y∂/∂y) exp(z∂/∂z) would lead to enormous technical problems. A better alternative consists in us...
متن کاملApplication of Laguerre Polynomials for Solving Infinite Boundary Integro-Differential Equations
In this study, an efficient method is presented for solving infinite boundary integro-differential equations (IBI-DE) of the second kind with degenerate kernel in terms of Laguerre polynomials. Properties of these polynomials and operational matrix of integration are first presented. These properties are then used to transform the integral equation to a matrix equation which corresponds t...
متن کاملThe Operational matrices with respect to generalized Laguerre polynomials and their applications in solving linear dierential equations with variable coecients
In this paper, a new and ecient approach based on operational matrices with respect to the gener-alized Laguerre polynomials for numerical approximation of the linear ordinary dierential equations(ODEs) with variable coecients is introduced. Explicit formulae which express the generalized La-guerre expansion coecients for the moments of the derivatives of any dierentiable function in termsof th...
متن کاملConstruction and implementation of asymptotic expansions for Laguerre-type orthogonal polynomials
Laguerre and Laguerre-type polynomials are orthogonal polynomials on the interval [0,∞) with respect to a weight function of the form w(x) = xαe−Q(x), Q(x) = m ∑ k=0 qkx , α > −1, qm > 0. The classical Laguerre polynomials correspond to Q(x) = x. The computation of higher-order terms of the asymptotic expansions of these polynomials for large degree becomes quite complicated, and a full descrip...
متن کامل